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An anisotropic power spectrum will have a clear signature in the 21cm radiation from high- 
redshift hydrogen. We calculate the expected power spectrum of the intensity fluctuations in neutral 
hydrogen from before the epoch of reionization, and predict the accuracy to which future experiments 
could constrain a quadrupole anisotropy in the power spectrum. We find that the Square Kilometer 
Array will have marginal detection abilities for this signal at z ~ 17 if the process of reionization 
has not yet started; reionization could enhance the detectability substantially. Pushing to higher 
redshifts and higher sensitivity will allow highly precise (percent level) measurements of anisotropy. 

I. INTRODUCTION 

Modern cosmology assumes that the universe is both homogeneous and isotropic when averaged on scales larger than 
100 Mpc. Observations, in particular that of the cosmic microwave background (CMB), have established homogeneity 
and isotropy as valid assumptions if we interpret them statistically. Put another way, "the density of matter may differ 
from one point in the Universe to another, but the distribution of matter is described as a realization of a random field 
with a variance that is everywhere the same and the same in every direction" [T]. With the accumulating precision 
in the measurement of the CMB temperature field and the soon expected polarization data from the Planck satellite, 
these assumptions can be tested. Ackerman, Carroll and Wise [2] have considered the possibility that rotational 
invariance was broken in a scale invariant way during inflation, by an effect that has since disappeared. Their paper 
suggests a possible explanation for hints of anomalies in recent CMB measurements, many of which are carefully 
discussed, and dismissed as statistical fluctuations in a recent analysis by the WMAP team [3] . Such a preferred 
direction would be an important clue to the physics of the early universe. PuUen and Kamionkowski [1] have developed 
CMB power spectrum statistics to detect a direction dependence in the fluctuations in temperature or polarization. 
Groeneboom et al. [3] have revisited the model of Ackerman et al. to include polarization and various systematic 
effects. 

While the CMB has been a tremendous resource for exploring the physics of the early universe, we are approaching 
the limits of the information that is encoded: the temperature fluctuations measurements expected from the Planck 
satellite [5] will be primarily limited by cosmic variance (limits on our ability to characterize variances due to the 
finite number of independent and informative measurements we can make) in the entire regime where the primary 
CMB fluctuations are largely unconfused by astrophysical foregrounds. Comparable measurements in polarization 
would not add a tremendous amount of new information, with the important exception of a possible detection of the 
signature of gravitational radiation generated in the early universe. 

Measurements of large scale structure can in principle provide much more information about the potential fluctua- 
tions; the CMB is largely a two dimensional surface, while surveys can probe the three dimensional structure of the 
universe. PuUen and Hirata [51 have used luminous red galaxies in the Sloan Digital Sky Survey to put the strongest 
limits on any possible anisotropy, constraining any possible quadrupole anisotropy in the power spectrum to be less 
than 40%. 

On large scales, the fluctuations are small enough such that they are expected to evolve according to linear theory, 
but small scales are highly non-linear and difficult to model. The scale at which this transition occurs defines the 
smallest scale which can be readily used for measuring cosmological parameters, which places limits on the number 
of independent modes of the fluctuations which can be measured. As structure grows in the universe, this scale 
progressively moves to larger scales. To explore a large number of modes, it is therefore helpful to make measurements 
at high redshift. 

A promising technique for measuring the three-dimensional structure of the universe at high redshift is to measure 
the 21cm hyperfine transition of neutral hydrogen. In this paper, we investigate the signature of a preferred axis for 
the fluctuation power in this redshifted 21cm emission. 

Neutral hydrogen is the dominant form of baryonic matter in the "dark ages" , i.e. before star formation, hence we 
expect redshifted 21cm radiation to reach us from all directions in the sky. The intensity of this radiation tells us 
about the distribution of neutral hydrogen in the universe, as a function of both angular coordinates in the sky and 
redshift. Thus, in contrast to the CMB, 21cm surveys can probe the three-dimensional distribution of matter in the 
universe (see [7] for an in-depth review). 
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As a benchmark experiment, we assume that it wiU be possible to make arcminute scale measurements with mK 
sensitivity in a bandwidth of order 1 MHz. This is a somewhat optimistic forecast, but not far beyond what is 
projected for the Square Kilometer Array [H]- It is expected that ~ 250 square degrees could be surveyed down to 
frequencies of 70 MHz to a sensitivity of 10 mK, with thousands of frequency resolution elements in a band. At these 
frequency and spatial resolutions, the smallest accessible scales will be of order a few comoving Mpc. We also consider 
the benefits of a next generation experiment, an idealized version of an Omniscope (18j . that can operate at slightly 
lower frequencies (to push to higher redshifts), but has a wide field of view. This is a challenging experiment, as the 
noise rises steeply with lower frequency, but requires only a modest expansion in frequency coverage and a feasible 
increase in sensitivity. 

In what follows below we focus on the signal from before reionization. This simplifies the analysis and is robust 
to the details of how the first stars formed, but is neglecting a possibly large signal The 21cm spin temperature 
during reionization could be substantially different from the CMB temperature in the presence of a UV background, 
providing a substantial amount of statistical power jTUl, especially at the lower range of redshifts we consider, z ~ 17. 
The details of reionization are sufficiently uncertain at this point that it is difficult to predict the utility of this epoch 
for precise measurements of the matter power spectrum. It is possible that coupling to Ly-a photons has already 
coupled the 21-cm spin temperature to the gas temperature at z ~ 17, but there are currently no observational 
constraints on this epoch, and there is no exhaustive theoretical search of parameter space that has determined the 
most likely physical conditions at z 17. In the sense of signal-to-noise, this would be a large boost over what is 
assumed in what follows, but it comes at the cost of increased astrophysical uncertainty. 



II. THE 21 CM BRIGHTNESS TEMPERATURE AND ITS POWER SPECTRUM 



We can consider the equation of radiative transfer along the line of sight through a hydrogen gas cloud and then 
compare the intensity of the 21 cm coming through the cloud to a hypothetical "clear view" of the 21 cm radiation 
from the CMB radiation. This difference in intensity is expressed through the brightness temperature [7] : 
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where we have taken the values ^lo = 2.85 x lO^^^^ s"\ Hq = 73 km s"^ Mpc"\ 1^21 = 1420 MHz, = 0.0425, 
n,n = 0.26. 

The brightness temperature in eq. [ijis a function of redshift for a chosen line of sight. Since redshift can be related 
to the radial distance, we can think of the brightness temperature as a function of x. Redshift z is related to the 
comoving radial coordinate x through 
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For 1 < z < 3000 we have a matter dominated universe with H{a) = Ho^yfl„l/a^ . Here and throughout, we assume 
that the peculiar velocities are non-relativistic. 

After recombination Compton heating through the residual free electrons keep the kinetic cosmic gas temperature 
Tx and the spin temperature Tg, equal to the CMB temperature T^. As illustrated in ref. [7] fig. 6, the cosmic 
gas begins to decouple from the CMB at around z ~ 300. At this point the cosmic gas is dense enough that the 
spin temperature Ts is cooled from the CMB temperature to the kinetic cosmic gas temperature Tk by collisions. 
By z ~ 100 the cosmic gas begins cooling adiabatically as Tk = 0.02 K(l -I- z)^. As the cosmic gas expands 
collisions between hydrogen atoms become less efficient at cooling the spin temperature. Photon interactions then 
slow the cooling of the spin temperature, an it is eventually again driven towards the CMB temperature. By the time 
reionization begins at around redshifts of 20, Ts ^ T^, making the 21cm signal relatively small. In this our analysis 
we are particularly interested in redshifts z, with 15 < z < 35 when T5 < T-y. 

We now consider fluctuations in the 21 cm brightness temperature by defining an average of dTb over all angular 
coordinates and defining S2i{x) through 



sn{x) = 5n{z){i + S2i{x)) 



(4) 
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Fluctuation in the brightness temperature arise through fluctuations in the baryon density Sf,, the neutral fraction 
6xHn fluctuations in the spin temperature Sts^ and the line of sight peculiar velocity gradient (5g„. 

We are interested in z much after recombination but before reionization. During this time the neutral fraction 
fluctuations are so small that we can just take xhi — 1- 

The fluctuations in spin temperature can be related to fluctuations in the CMB temperature St^, the cosmic gas 
kinetic temperature St^ and the baryonic density fluctuations St- This is because the spin temperature is determined 
solely by the temperatures and as long as UV scattering is negligible, which is true before reionization. The 
relationship between spin and kinetic gas temperatures is expressed through the collision coefficients Xc which describe 
the rate of scattering among hydrogen atoms and electrons [TTl [T^ : 

^_7^\ ^^^^ 
Ts J l + Xc \ Tk J 

where /csT^ = ?i27ri/2i = 0.0682 K. Fluctuation in the CMB temperature are so small that they can be ignored. 
Combining the remaining contributions to S21 we calculate that at linear order 

S2i{x) = l3b{z)Sb{x) + I3tk{z)Stk{x) -6oy{x), (6) 

The coefficients f3i, and /3tk are 

^^"^^ Tk-T,^ x,{l + x,)\'^ dlnTK^ ' dlnTK)' ^' 

Now is related to Sf, via a proportionality constant that depends only on the redshift distance I13|: 5tk{x) = 
g{z)5b{x). Fig. 2 in ref. [I3j gives the proportionafity constant g{z) for redshift z between 10 and 1000. We define 

p{z) = (ib[z)+g(z)(3TAz) (9) 

and hence 821 in eq. |6] becomes 

52i{x) = p{z)5b{S) - 59y{x) (10) 

For 15 < z < 35, P{z) is approximately 1.6 to within 3.5% accuracy. 
The line of sight peculiar velocity gradient 



H{z)/{l + z) c 



(11) 



introduce redshift space distortions. For scales large enough for linear theory to hold. Kaiser [T3] has shown that 
Sdv{k) — —f{k ■ rsY5matter{k) . Lahav et al. [15] explain that /(z) — [r2m(l + z^ {Hq/ H{z))^]^-^ which is essentially 
unity for the matter dominated universe at the redshifts we are considering. Furthermore we assume that 5b follows 
^matter for the scales we are considering so that 

5dy{k)^-{k-fs?k{k) (12) 

We define the power spectra via 

(^2l(fc)^2l(fc')> = {2TTfS^^\k + k')P2l{k), (13) 

{Sbik)Sb{k')) = {27:fS^-^\k + k')Pbik), (14) 

By writing a spherical harmonic expansion of the 21 cm fluctuations, S2i{x) = 52i{Xi 2;) — X); m ^im{x)Yim{S:) we can 
construct the angular power spectrum 



j^3iikx)3i'{kx')YL{k)Y,„Ak)P2i{k) (15) 
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III. THE ANGULAR POWER SPECTRUM IN THE FLAT SKY DISTANT OBSERVER 

APPROXIMATION 

We now calculate the 21 cm angular power spectrum in the flat sky approximation distant observer approximation. 
We look at a tile of sky at the north pole (in the z direction), at an average distance xoj thickness Ax, and angular 
size Ae x Ae. The tile is not too thick, so that the constants in redshift z, such as g(z) and fiiz) do not vary much 
and we can treat them as constants. In particular from now on we drop the z dependence and write simply /3. As 
discussed in the previous section, for the redshifts of interest in this paper /? « 1.6. We mask the &{x) so that it 
is non zero only on that tile. This permits us to extend our integrals beyond the tile's dimensions to simplify their 
evaluation. Incorporating Kaiser's line of sight approximation for red shift space distortions allows us to write 

^2i(fc) = (/3+(fc-r,)2) ^b(fc) (16) 

where we take the line of sight fs — z. 

We follow [16] for the definition of the a's, and applying it to our case, for large I we define a 2-D vector I — 
(I cos 4)1 , 1 sin (pi ) , such that 

~ J d\e<^~''-^S2iix,i) (18) 

so that 

X X X X J 27r V H + iyx^J X 

In order to get rid of the 1/x^ dependence which prevents us from having a diagonalized correlation matrix in 
momentum space we will define b = x^o,{hx) ^^id consider 6 as a function of k = l/x and x instead of I and x- 
We then Fourier transform in x and define a 3-D momentum vector q = (R, fca) 

mh\<i)) = {27Tf6^^\q + cf)P2i{q) (20) 
= {27rf S'i^\q + ^){p + {q-zf)'p,{q) (21) 

In doing the scaling by distance to make the signal covariance matrix diagonal, the noise covariance will become 
more complicated. However, the noise covariance matrix is expected to already be complicated due to non-trivial 
foregrounds. 

IV. FORECASTING THE ANISOTROPIES 

We follow PuUen-Kamionkowski [T] and expand the anisotropic Pb{q) in terms of spherical harmonics 

PM = Viq)[l + E 9LM{q)YLMm (22) 

To estimate the constraints possible on the ql.m^s, we calculate the Fisher matrix. As in jjj we will take the gL,M{q) 
to be constant in q. We will first do the case without noise, and then we will consider the case with noise. 
The Fisher matrix is 

F = ^Tr[C„,,,, C-i(C,3^ )tC-i] (23) 



We will consider a quadrupole anisotropy, i.e. L=2. Eq. 21 is the correlation matrix C^^ and 



C,,,,,, = {2n)H'^j^\q + q') (/3 + {q ■ zffv{q)YLMq) (24) 
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so that 

C,,.^.C-\C.,^^^)^C-^ = {2.)H^^[^+^)—^^^^^^^^^ (25) 

[1 + l^M" 9L,M"yL,AI"{q)r 

Taking the trace of the above means setting q' = —q and integrating with measure j^^- This leads to {2Tr)^S^{0) 
which equals the volume of our thick tile. Thus we have 

Fm,M' = 7,{Vol) / ——^j—^ — (26) 

The measure / d^q = q^dq J dilq can be regulated with the knowledge that our volume and resolution are both 
finite; we approximate the volume as a pixelized lattice with a total of = Nt pixels. Thus 



oo 



Thus we have 



(27r)3 7o " " {2tt)3 L L J \L J 3 3 Vol 



Nt f ,^ >L.M(Q)i2,M'(9) 



Fm^m' - dQq ' '— — -- (28) 

6 J [^ + T,M" 9l.M"Yl.M"{(1)\ 

Since we expect the gL.M" to be small we evaluate the integral with gL,M" = 0. By the orthonormality of the spherical 
harmonics we get the Fisher matrix without noise: 

Nt 

[-FAf,Af'] noiseless = —^5m,M' (29) 

This is trivial to invert. 

Now we include the beam size and instrument noise. The correlation matrix given in Eq. |21| becomes 

C,-,,-, = {2^f5'p{q + ,^)P2M (30) 
= {2nf5^;^\q + q') [(/? + {q ■ zff P,{q) exp{-aW) + Pn] (31) 

This differs from the noiseless case by the addition of two new parameters, ctb and P„. We discuss the meaning of 
each in turn. 

We assume a Gaussian beam characterized by a beam width as- This beam width is related to both the radial 
direction (i.e. frequency resolution as given by Eq.js]) and the transverse direction angular resolution. While in general 
these two resolutions, for the radial and ctt for the transverse, can be different, an experiment well constructed 
for detecting an anisotropic matter spectrum would choose them so that a^P + a'^k^ = o"^?^- 

We model the instrument noise as white, uncorrelated between pixels and uniform in momentum with value P„. 
We relate the noise per pixel temperature T„ to P„ through 

(TJ Snf = J exp(-a|qV(4^)) = ^ (32) 

Here we have incorporated that cutoff as a Gaussian characterized by the beam width, but any reasonable cutoff in 
the momentum volume would give similar results. Varying the power-law index of the noise (to 1/fc, for example), 
affects the forecasts. Such considerations can straightforwardly be included in our analysis and we discuss them in 
the conclusions. 

The derivative of the correlation matrix with respect to the 5L=2,Af's is 

C,g,,,, = {27rf6^^\q + q') (/3 + (g • z)^)' V{q) exp{~aW)yL,Mi9) (33) 
and the Fisher matrix becomes 

'M,M' = ni^ol)-—^ / d^q- ^ (34) 



2' '(27r)3 



1 + Ea/" 9L,M"YL,M"{q) + (/3+(q.£)2)^-p(^)cxp(-^2g2) 
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FIG. 1: Matter power spectra for z = 17.5 (thick dashed) and z = 30 (thick solid). The lower z curve is higher simply due to 
the growth of structure in linear theory. The thin blue curves show the noise per In interval in k assuming the experimental 
parameters outlined in the text. These curves include the impact of both experimental noise and the redshift evolution in the 
mapping between density fluctuations and temperature fluctuations. Substantial star formation at high z could lower the noise 
curve at 2 ~ 17.5 by up to 4 orders of magnitude. 



Again, since we expect the gL,M" to be small we evaluate the integral with gh.M" ~ 0. We are particularly interested 
in forecasting for the smallest redshifts before ionization, i.e. 15 < z < 35 which is where we can take /3 « 1.6. For 



these redshifts (/3 + (q • z)^ 
by 1/2 so that (/3 + (g • zf 



will be between 2.56 and 6.76. We approximate our Fisher matrix by replacing {q ■ 
is replaced by 4.41. 
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(35) 
(36) 
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We evaluate numerically the integral for two cases. Case I will correspond to the Square Kilometer Array target 
values. Case II will correspond to the Fast Fourier transform telescope (FFTT) described in [TS] . 

From the thermal noise per visibility as given by Morales [19j in eq. 7, we arrive at the thermal noise per pixel for 
the brightness temperature given by 



^sys ^diffractii 



-'desired 



(38) 



where Tgyg is the system temperature, B is the bandwidth, and r is the total observing time, ^diffraction is the diffraction 
limited resolution X2i{l + z)/\/Ae (A^ is the effective antenna area) and ^desired is the 1 arcminute resolution we desire. 
The angular resolution is assumed to be tuned by a dilution of the array from being fully compact by a simple scaling 
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of all baseline positions by a fixed amount. The system temperature is given by ARCADE 2 '50' as 



Tsys = 1.26 K (1 + z) 



1 GHz\ 

1^21 J 



2.6 



(39) 



Putting this together we arrive at an expression for the noise as a function of rcdshift assuming that an array is simply 
scaled to maintain constant angular resolution: 



12 niK 



1 + Z\3.85r ^ 



V(r/104hr)(^,/kr 



1 



1/2 



(40) 



The noise is a steeply rising function of redshift, but the signal is also a strong function of redshift in the dark ages 
just before reionization. For example, if star formation has not started in earnest by z ~ 17, we can see from the 
table below that the mean 21cm brightness temperature at z = 30 is 30 times larger, while the noise is higher by only 
a factor of 8. The growth of fluctuations boosts the lower z signal by a factor of 1.7, but that still leaves more than 
a factor of 2 higher sensitivity at z ~ 30 compared to z ^ 17. 

For both our cases we will consider 10 000 hours of total observing time. For case I, the SKA, we consider a line 
of sight depth of redshift z = 15 to 20, i.e. 374 Mpc/h. The effective antenna area is 1 sq km. This gives an average 
noise per pixel of about T„ = 8 mK. The total area of sky observed is taken to be 250 degrees. 

For case II, the FFTT (or "Omniscope"), we consider a line of sight depth of z = 25 to 35, i.e 346 Mpc/h. The 
effective antenna area is 100 sq km, the limit at which earth curvature could be a problem [18] . This gives an average 
noise per pixel of about = 6 mK. The total area of sky observed is taken to be 3000 degrees, to account for edge 
effects not allowing a full 1/2 sky. 

These two cases correspond to a central redshift values of z = 17.5 and z = 30 respectively. Because of our finite 
box size we impose a lower limit in the dq integration corresponding to 27r divided by the smallest dimension of the 
box, which is the line of sight direction. We have for the z G [15,20] and z G [25,30] a lower momentum cutoff 
of 0.0168 h/Mpc and 0.0181 /i/Mpc, respectively. In the numerical evaluation of the integral we have checked that 
replacing this lower momentum cutoff by zero does not significantly affect the result. 

We take our beam width to correspond to the length scale of 1 arcminute at the comoving distance of interest. 
For z—17.5 and z=30, the comoving distance is 9020 Mpc/h and 9650 Mpc/h, respectively, and the beam width is 
2.63 Mpc/h and 2.81 Mpc/h, respectively. 

We use the linear matter power spectrum as provided in the LAMBDA CAMB Web Interface Toolbox ^7j. The 
power spectra can be seen in Figure [T] for the different redshift choices, along with the noise curves for the assumed 
experimental parameters. For the high redshift case, the signal to noise is much higher for two reasons: the assumed 
noise level is slightly smaller in power (the increased noise at higher z is assumed to be more than offset by a larger 
collecting area), and the expected cosmological mean signal (in mK) is substantially larger, leading to a significantly 
stronger constraint on the matter power spectrum for the nominal higher redshift experiment. 

We summarize our results for the two cases in the following table. 



Case 


z 


sn 

[mK] 


r„ /pixel 
[mK] 


beam width 
[Mpc/h] 


angular size 
[sq. deg.] 


vol 

[{Gpc/hf] 


I 

[{h/Mpcf] 


T[g2M] = F-^/-' 


I 


15-20 


-0.150 


8 


2.63 


250 


2.32 


4.11 X 10-" 


4.6 


II 


25-35 


-4.42 


6 


2.81 


3000 


29.45 


2.04 X 10-6 


0.006 



V. CONCLUSIONS AND DISCUSSION 



The results in the above table mean that 21 cm surveys should be able to constrain the values of (72M to about 
l/VO.0476 = 4.6 for case I and l/\/300()0 = 0.006 for case II. If instead of white noise, we allow for a Pn{q) ^ 1/?, 
then the constraints on the values of g2M are 43 and 0.03 for cases I and case II, respectively. 

While the noise per mode is quite high, the large volume allows a huge number of modes to be measured. For 
example, 1 (Gpc/h)^ of volume provides (in principle) more than 3 x 10^ measurements with ~ Mpc/h resolution. 
Even with such a large number of samples, we see that this is a challenging measurement with planned upcoming 
experiments such as SKA; to accurately measure this signal with 21cm experiments will necessitate deeper maps, to 
push the signal to noise per mode high enough to reach the sample variance limit. 

One way to increase the signal to noise is to measure the signal at higher redshift, where the pre- reionization 
hydrogen spin temperature difference from the CMB is larger. At z ~ 28 the mean temperature difference is 30x 
larger than at z ~ 17, in the absence of reionization effects, while the sky noise is higher by less than a factor of 10. 
It has been suggested that reionization could boost the z ^ 17 signal by up to two orders of magnitude Pj if the 
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reionization process doesn't compromise the signal of interest, this would allow SKA to have a measurement of the 
anisotropy to better than ten percent. 

The promise of future 21cm experiments is evident from the large signal to noise that is possible from a sufficiently 
sensitive experiment. This is in contrast to CMB experiments, where cosmological information is now largely limited 
by the finite number of modes that can be measured within the boundaries imposed by causality. For future 21cm 
experiments, limits will be set by experimental capabilities, but precise measurements at high rcdshift (i.e., lower 
frequencies) will allow powerful constraints on fundamental physics. 
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